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ABSTRACT 

The twin peaks in the nucleus of M31 have been interpreted by Tremaine as a 
thick, eccentric, disk of stars orbiting a massive dark object; the required align- 
ment of the apoapsides of the stellar orbits could be maintained by self-gravity, 
and the whole structure might be a discrete, nonlinear eigenmode. The pattern 
speed of this mode could, in principle, be determined by the Tremaine- Weinberg 
(TW) method, which requires measurements of the surface brightness, and ra- 
dial velocity along a strip parallel to the line of nodes. However, spectroscopic 
observations along the line of nodes are not available. We propose a variant of 
the TW method, which exploits a basic feature of the eccentric disk model, to 
extract estimates of the pattern speed from Hubble Space Telescope spectroscopic 
data, taken along the line joining the two peaks. Within limitations imposed by 
the data, we estimate that the pattern rotates in a prograde manner and, for 
an assumed disk inclination of 77°, the pattern speed \VL P \ < 30 kms _1 pc _1 , or 
period more than 200, 000 years. 

Subject headings: galaxies: individual (M31) — galaxies: kinematics and 
dynamics — galaxies: nuclei 

1. Introduction 

The nucleus of M31 was first resolved by the Stratoscope II balloon-borne telescope 
(Light, Danielson, and Schwarzschild 1974), which showed that the peak brightness was 
displaced relative to the center, as inferred from the isophotes of the outer parts of the 
galaxy. This was confirmed, and extended by Hubble Space Telescope (HST) observations, 
which revealed two peaks in the brightness, separated by 0".49 (Lauer et al. 1993, King, 
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Stanford, and Crane 1995). Ground-based, as well as HST spectroscopy, have probed the 
structure of the radial velocities and velocity dispersions, in increasing detail, along many 
strips across the nuclear region (Dressier and Richstone 1988, Kormendy 1988, Bacon 
et al. 1994, van der Marel et al. 1994, Gerssen, Kuijken, and Merrifield 1995, Statler et 
al. 1999, Kormendy and Bender 1999). These provide evidence for the presence of a massive 
dark object (MDO), which could be a supermassive black hole, of mass ~ 3 x 10 7 M , 
located very close to the fainter peak (P2). The dynamical center of the nucleus is believed 
to coincide with the center of the isophotes of the bulge of M31; this point has been 
estimated to lie between the two peaks. 

Tremaine (1995) proposed that the nucleus could be a thick eccentric disk, composed 
of stars on nearly Keplerian orbits around the MDO, with their apoapsides aligned in the 
direction toward the brighter peak (PI); the brightness of PI is then explained as the 
increased concentration of stars, resulting from their slow speeds near their apoapsides. He 
also showed that this model is broadly consistent with the kinematics, as inferred from the 
spectroscopic observations of Kormendy (1988), and Bacon et al. (1994). Recent work has 
not only produced further support for Tremaine's model (Statler et al. 1999, Kormendy and 
Bender 1999), but has stimulated variations on the basic theme (Statler 1999). Tremaine 
also suggested that the alignment could be maintained by the self-gravity of the disk, 
wherein the eccentric distortion could arise as a discrete, nonlinear eigenmode, with some 
nonzero pattern speed (fl p ), equal to the common apsidal precession rate. The dynamical 
question is yet to be resolved in a self-consistent manner, although explorations of orbits in 
model potentials have identified a family of resonant, aligned loop orbits, which could serve 
as building blocks (Sridhar and Touma 1999, Statler 1999); reasonably faithful reproduction 
of the nuclear rotation curve adds some degree of confidence (Statler 1999). If the nuclear 
disk is indeed a steadily rotating, nonlinear eigenmode, what is fl p ? 
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Tremaine and Weinberg (1984, hereafter TW) invented a method of estimating the 
pattern speed of a barred disk galaxy, that uses measurements of the surface brightness, and 
radial velocity along a strip parallel to the line of nodes (defined as the line of intersection 
of the disk and sky planes). That this methods works was proved when the pattern speed 
of the bar in NGC 936 was estimated by Merrifield and Kuijken (1995). Unfortunately, the 
radial velocity measurements of the nucleus of M31 (see references above on spectroscopy) 
are available on strips that, either do not coincide with the line of nodes, or possess too 
poor an angular resolution for direct application of the TW method. In this Letter, we 
show that the HST observations of Statler et al. (1999, hereafter SKCJ), together with 
Tremaine's (1995) model, can be used to estimate Q p . 

2. A modified Tremaine— Weinberg method applicable to M31 

We briefly recall TWs derivation of their kinematic method. Let us assume that 
the disk is razor-thin, with a well-defined pattern speed, Q p . The plane of the disk is 
assumed to be inclined at angle i with respect to the plane of the sky. Let (x,y), and 
(r, 0) be cartesian and polar coordinates, respectively, in the plane of the disk, with the 
origin coinciding with the center of mass of the system (disk plus MDO). Let the cartesian 
coordinates in the sky plane be (X,Y) = (x, ycosi) ; the x and X axis coincide with the 
line of nodes. The disk is assumed to rotate steadily, hence the surface brightness of stars, 
E(x, y,t) = E(r, <fi — Q p t) . The surface brightness is also assumed to obey a continuity 
equation, without a source term. Let v be the (mean) velocity field in the disk plane. The 
continuity equation can be manipulated to yield, 

V-[(v-u)E] = 0, (1) 

where u = Q p (—yx + xy) . TW proceed by integrating equation (1) over x from — oo to 
+oo. Assuming E — > sufficiently rapidly as \x\ — > oo, and integrating over y from y to 
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oo yields Qpf^dxT.x = J^dxHv y . Noting that the sky brightness, S s = E/ cos?, 
and the radial velocity, V\\ = v y sini, the integrals may be expressed in terms of observable 
quantities; hence fl p can be estimated when V\\ has been measured on a strip parallel to the 
line of nodes. 

SKCJ measured V\\ along the P1-P2 line, which is inclined by about 4°, in the 
sky plane, to the line of nodes. Therefore V\\ is available on a strip that makes an 
angle, 9 ~ 4°/ cos 77° ~ 18°, with the x-axis, in the disk plane; the TW procedure 
needs some modification, before Q p can be extracted. Let x' — x cos 9 + y sin 9 and 
y' = —x sin 9 + y cos 9 be the rotated cartesian coordinates. Let us also denote the surface 
brightness by T,'(x',y',t) . The SKCJ measurements are along the strip y' — , that passes 
through the origin. Equation (1) expresses a relation that is invariant under rotation 
of cartesian coordinates. Hence application of the TW procedure provides an identical 
relationship between the integrals over the strip, defined by y' = : 

/oo poo 
dx'E'x' = / dx'Z'v' y (2) 
-oo J — oo 

We now express the integrals in terms of observable quantities: x' = X cos 9 + Y sin9/ cosi 
and = y' = —X sin 9 + Y cos 9/ cosi can be used to eliminate Y, giving x' = X/ cos 9 . As 
before, Y>'(x',y',t) = cosi S S (X, Y, t) . However, the radial velocity, 

V\\ = sin i (y' y cos 9 + v' x sin 6>) , (3) 

depends on v' y as well as v' x . Hence it is, in general, not possible to express v' y in equation (2) 
in terms of V\\ alone (the TW method finesses this problem because 9 = kills the 
contribution to V\\ from the v' x term). However, we are able to make progress by recalling 
the basic elements of Tremaine's (1995) model, //the nuclear disk is largely composed of 
nearly Keplerian ellipses, with their apoapsides aligned along the P1-P2 line, and we choose 
9 such that our strip lies along this line, then the strip intersects all these ellipses at right 
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angles. Thus v' x — , and we recover the familiar TW relation, 

/oo /-OO 
dXY> s X = I dXT. a V\\. (4) 
-oo J — oo 

It should be noted that the value of 9 drops out of equation (4). Moreover, v' x = even 
when the apsides of the ellipses precess. Below we use HST photometry and kinematics to 
estimate Q„ . 



3. Q p from HST photometry and spectroscopy 

HST photometric data, reported in Lauer et al. (1993), was kindly supplied to us by 
Prof. Ivan King. Figure 1 shows the (sky) surface brightness along the P1-P2 line. The 
contribution from the bulge of M31 was subtracted using two different fits to the bulge 
brightness, namely a Nuker fit as described in Tremaine (1995), and a Sersic fit as described 
in Kormendy and Bender (1999). The disks so obtained will henceforth be referred to as a 
"T-disk", and a "KB-disk", respectively. 

SKCJ observed the stellar kinematics along the P1-P2 line, using the //48 long-slit 
spectrograph of the HST Faint Object Camera. We obtained radial velocities along the 
P1-P2 line, including the errors on them, from their "de-zoomed" rotation curve (given 
in Figure 3, as well as Table 1 of SKCJ); these are displayed in Figure 2a. The integrals 
in equation (4) need to be computed, with upper and lower limits symmetrically displaced 
about the center of mass of the disk plus MDO. Kormendy and Bender (1999) have 
determined the center of mass to be displaced by 0".098 from P2 toward PI, and we use this 
value in the computations for Figures 2b and 2c. The stated errors on the radial velocities 
were used by us to generate 300 random realizations. For each of these realizations of the 
rotation curve, we evaluate the integrals for five different limits, ranging from ±0".8 to 
±1".2 . As is clear from equation (4), only the combination, Q p sin 2 can be determined. 
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It is a common assumption that the nuclear disk of M31 has the same inclination, to the 
sky plane, as the much larger galactic disk of M31, which is inclined at ~ 77° . We wish to 
state our results independent of this assumption, so we present estimates of the quantity, 
£l p = (smi/ sin 77°) fi p , in Figures 2b— 2d. 

Figures 2b and 2c display the estimates of Cl p so obtained, together with la error bars, 
as a function of the limits of the integrals, for the T-disk and KB-disk, respectively. It is 
evident that the estimates of pattern speed do not vary much when the limits of integration 
lie between ±1".0 and ±1".2 . We also explore the dependence on the position of the center 
of mass of the system, assumed in the computation of the integrals. Although we used the 
most recent determination, due to Kormendy and Bender (1999), it must be noted that 
earlier work (Lauer et. al. 1993, King, Stanford, and Crane 1995, Tremaine 1995) records 
smaller values, < 0".05, away from P2 toward PI. Figure 2d plots our estimates of Q p , 
for a range of values of the center of mass, with limits of integration fixed at ±1".0 . A 
negative value of Cl p means that the pattern is prograde; Figures 2b — 2d indicate that this 
is the most likely possibility, with the absolute value increasing with the separation between 
P2 and the center of mass of the system. This is reasonable, because a larger separation 
corresponds to a larger disk mass, which implies a greater contribution from the self-gravity 
of the disk, which is ultimately responsible for the (alignment and) precession of the disk. 
The errors on Q p are large, and a non rotating pattern cannot be ruled out with certainty. 
Below we quote representative bounds on the absolute value of the prograde pattern speed: 

sin 77° f 34 ± 8 km s" 1 pc _1 , T-disk; 
\%\<-^-\ (5) 
sm * I 20 ± 12 kms^ 1 pc~\ KB-disk. 

For each realization of the rotation curve, the zero-crossing point (the "rotation center" ) 

was determined by a third-order spline interpolation. We found the position of the rotation 

center to be displaced by 0".17 ± 0".01, from P2 toward PI. This should be compared with 

the value of 0".16 ± 0".05, quoted by SKCJ. We also tested for any systematic dependence 
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of Q p on the position of the rotation center, and found none. 

4. Conclusions 

Our estimate of the pattern speed of the nuclear disk of M31 should be qualified by a 
discussion of possible sources of errors, most of which are difficult to estimate quantitatively. 
SKC J calibrate velocities relative to an average over an 8 arcsec 2 aperture centered on the 
nucleus (Ho, Filippenko, and Sargent 1993), and assure us that the errors are likely to be 
small. More significant, perhaps, are the systematic errors in V», mentioned by SKC J; these 
are shown as open squares in Figure 2a. SKCJ used a slit of width 0".063, and this will 
introduce contributions to V\\ from nonzero values of v' x (see equation (3)). This effect is 
somewhat mitigated by cancellation between positive and negative values of v' x , and the 
fact that the width of 0".063 is of much smaller scale than the minor axis, projected onto 
the sky plane, of the smallest ellipse used by Tremaine (1995; see Figure 2a of his paper) 
to represent the nuclear disk. The limits of integration are necessarily finite in numerical 
computation, and we have resisted the temptation to include corrections by extrapolation 
of the S s and V\\ profiles. 

A basic assumption underlying the application of a TW-like method is that the 
surface brightness obeys a continuity equation, which would be valid for a stellar disk in 
the absence of star formation (or death). We expect the numbers of stars to be conserved, 
except possibly in the vicinity of P2, where the observed UV excess has been interpreted as 
contributions from early-type stars (King, Stanford and Crane 1995, Lauer et al. 1998). 
However, these stars do not contribute much to the photometric and kinematic data we 
have used. Tremaine (1995) argues that two-body relaxation is expected to thicken the 
disk, whereas we assumed that the disk was razor-thin. The original TW method is 
applicable to thick disks, so long as the streaming velocity normal to the disk plane is zero. 
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In addition to the assumption of zero normal streaming velocity, let us suppose that the 
three dimensional density, p, is symmetric about the midplane of the disk. The contribution 
to Vj| from v x arise from an integral along the (inclined) line of sight that runs through the 
thick disk. Consider two points along this line of sight that are equally displaced about 
the midplane of the disk. Aligned, nearly Keplerian orbits have flows such that p is equal, 
whereas v' x is equal and opposite at these two points; in this ideal picture, there is pair-wise 
perfect cancellation, and no net contribution to V\\ from v' x . In practice there should be 
some cancellation, because v' x will have opposite signs at two corresponding points, but 
there could be a net contamination from the unequal values of \v' x \ and p. 

Tremaine's original model, which was a reasonable fit to the then available photometry 
and kinematics, considered a non rotating disk, and it would be appropriate to inquire 
about the implications of a non zero pattern speed. A pattern that is prograde with angular 
speed, say, of 20 kms _1 pc _1 would contribute only about 35 kms -1 to the radial velocity 
at PI, which is about 250 kms -1 , according to the measurements of SKCJ. The maximum 
radial velocity quoted by Tremaine (1995) is less than 200 kms -1 , so a non zero pattern 
speed could still be accomodated. Our estimates do not rule out a non rotating disk, but 
we would like to offer a physical argument in support of a non zero pattern speed. For the 
disk plus MDO to be in a steady, non rotating state, the gravitational force on the MDO 
should necessarily vanish. Our (unpublished) numerical computations indicate that the 
force is indeed non zero. 

A limitation of our method is that it uses, in an essential manner, the assumption 
that most of the contribution to V\\ comes from orbits which intersect the measurement 
strip at right angles. In comparison, the original TW method does not rely on assumptions 
about the geometry of the mean flow; averaging over several strips, all parallel to the line 
of nodes, will improve estimates of the pattern speed, as Merrifield and Kuijken (1995) 
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demonstrated. Thus it is necessary to verify our estimates of Q p , by using the TW method 
on future observations of V\\, together with better photometry such as Lauer et al. (1998), 
along strips parallel to the line of nodes. An extremely useful set of observations that could 
be performed would be two-dimensional spectroscopy, similar to the work of Bacon et 
al. (1994), with the increased angular resolution that should be available in the near future. 

We are grateful to Prof. Ivan King for generously sharing with us the photometric 
data, an anonymous referee for pointing out a basic error in the original manuscript, and 
asking stimulating questions, and to R. Srianand for useful comments. NS thanks the 
Council of Scientific and Industrial Research, India, for financial support through grant 
2-21/95(II)/E.U.II. 
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Fig. 1. — V-band surface brightness along the P1-P2 line, with P2 located at X = 0, and 
PI at X = — 0".49. The solid line connects data points taken from Lauer et al. (1993). 
The dotted and short-dashed curves are the contributions from Nuker (Tremaine 1995) and 
Sersic (Kormendy and Bender 1999) bulges, respectively. The corresponding estimates for 
the surface brightness of the disk (T-disk and KB-disk) are given by the long-dashed, and 
dot-dashed curves. 



- 12 - 



200 - 



s 



o 

> 



-200 - 



3 -400 




-2 



20 



- 



o 



-20 



-40 
-60 



_L 



- 1 . 

X (arcsec) 



1 

(c) 



T 



KB-disk 



_L 



_L 



0.8 1 1.2 

disk radius (arcsec) 



o 




0.8 1 1.2 

disk radius (arcsec) 

(d) 



- 



o 



* -20 - 



-40 - 



T 



T 



T 



1 1 I 1 11 I 
□ T-disk " 
■ KB-disk - 



_L 



_L 



_L 



_L 



0.02 0.04 0.06 0.08 0.1 

center of mass (arcsec) 



Fig. 2. — (a) Radial velocity measurements taken from SKCJ. (b, c) Q p plotted versus the 
truncation radius of the disk, as measured from the center of mass of the system, for the 
T-disk and KB-disk, respectively; the center of mass was taken to be 0".098 away from P2 
toward PI. (d) fi p plotted versus the assumed location of the center of mass of the system, 
for a disk truncation radius equal to 1".0. 



-13- 
REFERENCES 

Bacon, R., Emsellem, E., Monnet, G., and Nieto, J.L. 1994, A&A, 281, 691 

Dressier, A., and Richstone, D. O., 1988, ApJ, 324, 701 

Gerssen, J., Kuijken, K., and Merrifield, M. R. 1995, MNRAS, 277, L21 

Ho, L. C, Filippenko, A. V., and Sargent, W. L. W. 1993, ApJ, 417, 63 

King, I. R., Stanford, S. A., and Crane, P. 1995, AJ, 109, 164 

Kormendy, J. 1988, ApJ, 325, 128 

Kormendy, J., and Bender, R., 1999, ApJ, 522, 772 

Lauer, T. R., et al. 1993, AJ, 106, 1436 

Lauer, T. R., et al. 1998, AJ, 116, 2263 

Light, E. S., Danielson, R. E., and Schwarzschild, M. 1974, ApJ, 194, 257 
Merrifield, M. R. and Kuijken, K. 1995, MNRAS, 274, 933 
Sridhar, S., and Touma, J. 1999, MNRAS, 303, 483 

Statler, T. S., King, I. R., Crane, P., and Jedrzejewski, R. 1. 1999, AJ, 117, 894 

Statler, T. S. 1999, ApJ, 524, L87 

Tremaine, S. and Weinberg, M. D. 1984, ApJ, 282, L5 

Tremaine, S. 1995, AJ, 110, 628 

van der Marel, R. P. et al. 1994, MNRAS, 268, 521 

This manuscript was prepared with the AAS IATfTjX macros v5.0. 



